Topologically distinct classes of valence bond solid states with their parent 

Hamiltonians 



Hong-Hao Tu and Guang-Ming ZhanjJ 
Department of Physics, Tsinghua University, Beijing 100084, China 

Tao Xiang 

Institute of Physics, Chinese Academy of Sciences, P.O.Box 603, Beijing 100190, China; 
Institute of Theoretical Physics, Chinese Academy of Sciences, P.O.Box 2735, Beijing 100190, China 



o 
o 

(N 



o 



> 
o 
in 
vn 
o 

^' 

o 

On 
O 



X 



Zheng-Xin Liu and Tai-Kai Ng 
Department of Physics, Hong Kong University of Science and Technology, Kawloon, Hong Kong, China 

(Dated: July 2, 2009) 

We introduce a general method to construct one-dimensional translationally invariant valence 
bond solid states with a built-in Lie group G and derive their matrix product representations. The 
general strategies to find their parent Hamiltonians are provided so that the valence bond solid states 
are their unique ground states. For quantum integer spin-S chains, we discuss two topologically 
distinct classes of valence bond solid states: One consists of two virtual SU{2) spin- J variables 
in each site and another is formed by using two SO{2S + 1) spinors. Among them, a new spin-1 
fermionic valence bond solid state, its parent Hamiltonian, and its properties are discussed in detail. 
Moreover, two types of valence bond solid states with 50(5) symmetry are further generalized and 
their respective properties are analyzed as well. 
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I. INTRODUCTION 

In recent years, the study of topological order has be- 
come a common issue in condensed matter physics and 
quantum information theoryjii^i^iii^i^ii Historically, this 
concept was proposed to describe fractional quantum 
Hall states,— which are incompressible quantum liquids 
with a finite energy gap to all bulk excitations. These 
new quantum phases of matter can not be described by a 
local order parameter with spontaneous symmetry break- 
ing. So the discovery of fractional quantum Hall effect 
brings great challenge to the Ginzburg-Landau theory, 
which is a corner-stone paradigm to characterize phases 
and phase transitions in condensed matter physics. From 
the viewpoint of quantum information theory, the ap- 
pearance of a long-range quantum entanglement plays an 
essential role in the topologically ordered states. How- 
ever, a general multipartite-entanglement measure that 
captures the most relevant physical properties is still 
lacking, because the number of parameters required to 
describe a quantum many-body state usually grows ex- 
ponentially with the particle number. 

The topological order appears not only in the two- 
dimensional fractional quantum Hall states but also 
in one-dimensional systems, for instance, the quan- 
tum integer-spin chains. In 1983, Haldane predicted 
that quantum integer-spin antiferromagnetic Heisen- 
berg chains have an exotic energy gap.^ Later, Affleck, 
Kennedy, Lieb, and Tasaki (AKLT)i^ found a family of 
exactly solvable integer-spin models with valence bond 
solid (VBS) ground states, and the presence of an ex- 
citation gap can be proved rigorously. In a spin coher- 
ent state representation, these VBS states share a strik- 



ing analogy to the fractional quantum Hall statesiii Al- 
though the two-point spin correlation decays exponen- 
tially, den Nijs and Rommelse^ found a hidden topo- 
logical order in the S ^ 1 VBS state, which is charac- 
terized by non-local string order parameters. For the 
5=1 VBS state, the long-range string order and the 
fourfold degeneracy in an open chain can be understood 
as natural consequences of a hidden Z2 x Z2 symme- 
try breaking ] ^^1^"^ For the standard integer-spin Heisen- 
berg models, the existence of spin-S'/2 edge states was 
also verified by quantum field-theory mappingsiS, and 
numerical calculations^—, which coincide with the VBS 
picture of the AKLT models. Experimentally, the VBS 
picture for S — I Haldane chain was supported by 
the electron spin resonance studieaiiii^ of the compound 
Ni(C2H8N2)2N02(C104) (NENP), the NMR imaginingi^ 
and the magnetic neutron scattering study2£ of the quasi 
one-dimensional material Y2BaNi05. 

The one-dimensional VBS states can be represented in 
a matrix-product form.— i^^iS^ Moreover, it was foundSi 
that density matrix renormalization group2^ (DMRG), 
the most powerful numerical method for one-dimensional 
quantum systems, converges to a matrix-product wave 
function as its fixed point. This important observation 
stimulates the formulation of the numerical techniques in 
one dimension by using the matrix-product variational 
wave functions. From a quantum information perspec- 
tive, the validity of the matrix-product variational ansatz 
depends on whether the true quantum ground states of 
the system obey an area law of entanglement entropy.— 
In this sense, the VBS states are only slightly entan- 
gled because their entanglement entropies have upper 
bounds even in the thermodynamic limit. Recently, these 



VBS states have received considerable attentions since 
they provide a playground to test the proposed measures 
of multipartite ent anglement . , 2 9 , 30 Xowards the poten- 
tial applications, it was suggested that the VBS states 
ensure measurement-based quantum computation;^ For 
the 5 = 1 VBS state of the AKLT model, Brenncn and 
Miyake'^^ have shown that a gap-protected measurement- 
based quantum computation can be performed within the 
degenerate ground states. 

In this paper, we will introduce one-dimensional trans- 
lationally invariant VBS states with a built-in Lie group 
G and present a general method to construct their par- 
ent Hamiltonians. For quantum integer spin-5' chains, 
we focus on two classes of VBS states. The local spin-S* 
states are formed by two virtual SU (2) spin- J variables 
in the first class and by two 5*0(25+1) spinors in the sec- 
ond one. To illustrate the method, we choose the 5" = 1 
fermionic VBS state with virtual spin J = 1 as an explicit 
example to find the parent Hamiltonian. Apart from the 
5 = 1 VBS state with virtual spin J = 1/2 and the 5 = 2 
VBS state with virtual spin J = 3/2 as the intersection 
elements, the VBS states of the two classes are shown 
to be topologically distinct to each other, which can be 
characterized by their edge spin representations in open 
chain systems. We also apply our method to investigate 
S'0(5) symmetric spin chains and discuss several VBS 
states with interesting properties. 

The outline of this paper is structured as follows. In 
Sec. II, we will introduce VBS states with a Lie group 
symmetry and derive their matrix product representa- 
tions. In Sec. Ill, we will focus on quantum integer- 
spin chains and study two topologically distinct classes 
of VBS states, including spin-S" VBS states formed by 
virtual SU{2) spin- J particles and by virtual 50(25-|- 1) 
spinor particles. Moreover, a spin-1 fermionic VBS state 
is extensively studied as an example and we construct its 
parent Hamiltonian explicitly. Sec. IV is devoted to the 
5*0 (5) symmetric VBS states and their physical proper- 
ties. In Sec. V, some conclusions are drawn. 



II. GENERAL CONSTRUCTION OF VBS 
STATES 

A. Matrix-product form 

We begin with a quantum one-dimensional chain with 
N lattice sites. In each site, the states (m = 

1, . . . , d) transform under a d-dimensional irreducible rep- 
resentation (IR) Gd of a Lie group G. Let us imagine that 
the physical Hilbert space is formed by two virtual iden- 
tical particles, whose internal quantum numbers {la)} 
{a = 1,...,D) transform under the _D-dimensional IR 
Gd of the same Lie group G. Here we require that both 
singlet representation G/ and IR Gd are included in the 
tensor product decomposition of two Go's. The first re- 
quirement means that Gd is a self- conjugate IR, i.e., the 
complex conjugate representation of Gd is equivalent to 
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FIG. 1: The schematic of the VBS states with a built-in Lie 
group G. Each dot denotes a virtual particle transforming 
under Gd irreducible representations of Lie group G. The 
solid lines represent valence-bond singlets formed by two vir- 
tual Gd irreducible representations on the neighboring sites, 
and the circles denote the projections of the virtual particles 
in each lattice site onto the physical Gd irreducible represen- 
tations. 

Gd- The latter requirement can be implemented by the 
projection operator onto the physical Hilbert spaced 

^-E E (1) 

m—l a, (3—1 

where is the Clebsch-Gordan coefficient defined by 

P^p = {Gd,m\GD,a;GD, P)- For VBS states in a peri- 
odic chain, each lattice site forms a valence-bond singlet 
I J) with its neighboring sites by pairing two virtual par- 
ticles (See Fig. [IJ. Thus, the wave functions of the VBS 
states can be expressed as 

I*) = (®f=li^fc,fc)|/)l2l^>23---|^)m> (2) 

where the valence-bond singlet |/) is given by 

D 
a,/3=l 

Here Qa,i3 = {Gi, I\Gd,Q:; Gd, P) is the Clebsch-Gordan 
coefficient of contracting two virtual Gd representations 
to form a singlet representation G/. 

In the present formalism, the VBS states can be easily 
written in a matrix-product form. Since P^"^ and Qa^p 
can be viewed as the matrix elements of D x. D matrices 
of and Q, the VBS states in Eq. ^ can be thus 
written as the following matrix-product form: 

1^-) = Tr(A[™iU[™^]---A[""l)|mi---mw), 

mi ■ ■■niN 

(4) 

where AI™! = pI^lQ is a L> x L> matrix. 

In periodic boundary conditions, the VBS states are 
invariant under lattice translation and transformation of 
the Lie group G by construction. Although no local order 
parameters can be found to characterize these states, the 
matrices can fully determine their physical proper- 
ties and render a "local" description. In open bound- 
ary conditions, edge states emerge at the two ends of 
the chain and then the matrix-product form of the VBS 
states is given by 

|«'a,/3)= J2 (^'''"^'^''"^'•■■^''""')a,/3|^l---^"iv), 
mi ---miv 

(5) 
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where the matrix indices a, (3 denote the edge states. 
These edge degrees of freedom are described by two frac- 
tionahzed particles transforming under Go representa- 
tion of the Lie group G. ActuaUy, the edge states and 
their representations are characteristic features of the 
VBS states, because they fully determine the local AI™! 
matrices. 

Another useful way to represent the VBS states is to 
use boson or fermion realization methods. To illustrate 
this method, we trace back to the tensor product decom- 
position of IRs of Lie algebras. According to the group 
theory, the physical states \m) under the exchange of the 
two identical virtual particles is either symmetric or an- 
tisymmetric, depending on Gd and Gd- Thus, the two 
virtual particles with fermion statistics create the anti- 
symmetric states, while the bosonic particles yield the 
symmetric ones. In some cases, there are still several 
channels with the same exchange symmetry and addi- 
tional projection has to be used to single out the phys- 
ical states in Gd- For example, the fermionic realiza- 
tion of a spin-2 VBS state with virtual spin J = 3/2 
was considered in Ref.'^^. In this case, both site-quintet 
states [S = 2) and site-singlet state {S = 0) are allowed 
for two spin-3/2 fermions on a single site and an extra 
projection can remove the unphysical site-singlet state. 
There also exists the Schwinger boson realization which 
symmetrizes several fundamental IRs to form higher di- 
mensional Gd's.— i^i^i^ In fact, all these boson/fermion 
realization methods play the role of (sometimes partially) 
deleting the unphysical states. 



B. Parent Hamiltonian: Locating the null space 

Following the spirit of the AKLT model, one can con- 
struct the parent Hamiltonians such that the VBS states 
in Eq. ^ are their unique ground states. It is most con- 
venient to work with the matrix-product form. For the 
matrix product states in Eq. ([4]), one can readily find 
that their reduced density matrix pi of I successive sites 
has a rank of at most. This suggests that the reduced 
density matrix pi of these VBS states always have null 
space for sufficient large I. These states are annihilated 
by the local projection operators supported in the null 
space. Hence, they are always exact zero-energy ground 
states of the translationally invariant Hamiltonians 

= (6) 

i 

where hi contains a sum of the positive semi-definite pro- 
jection operators supported in the null space from site i 
to i -\- 1 — 1. Previously, the parent Hamiltonians of the 
matrix product states for spin-ladder systems had been 
studied by similar methodsi^ 

Let us begin with the simplest cases with only nearest- 
neighbor interactions. Now the null space can be ob- 
tained from the VBS picture of these states. The Hilbert 



FIG. 2: The schematic of the "coarse-graining" process that 
converts the spins of successive sites to block spins. This 
procedure leads to a matrix product state with block spins, 
and the null space of a block-spin reduced density matrix can 
be identified. 

space of two neighboring sites can be divided into a di- 
rect sum of different IR channels according to the tensor 
product decomposition Gd® Gd- Once a valence-bond 
singlet of two virtual Gu's is formed, the remaining two 
particles of adjacent sites can transform under a direct 
sum of IRs resulting from Gd ® Gd- Therefore, the IR 
channels contained in Gd ® Gd but absent in Gd ® Gd 
constitute the null space in the 2-site reduced density 
matrix. 

The above steps to locate the null space can be em- 
bedded in a matrix-product formalism. Practically, we 
rewrite the matrix-product states in Eq. ([¥]) as 

|*> =Tr(gi52---gw), (7) 
where the local matrix gi is defined by 

g^^^Mj^^). (8) 

rrii 

To locate the null space, we resort to a "coarse-graining" 
procedure which converts the spins of adjacent sites to 
block spins. Since the VBS states are invariant under 
lattice translation, we can block the spins in sites 1 and 
2 as 

m-i ,m2 

= ^ B^o..Mg^c^^^Mg)^ (9) 

Gi2,Mg, 

where the D x D matrices ijI'^i^-^^Sl are given by 
5[Gi2,Afg]^ J2 A[™iU[™^l(Gi2,Mf2|Gd,mi;Gd,m2). 

mi ,m9 

(10) 

Here Gi2's distinguish the IRs of nearest-neighbor bond 
spin channels and |Gi2, are the states in IR channel 
Gi2. Correspondingly, {Gu, M^\Gd,mi;Gd,m2) is the 
Clebsch-Gordan coefficient to combine the states of G^'s 
into the states of G12. In the example of SU{2), G12 
denotes for the total bond spin S12 and —S12 < Mf2 < 
S12- After this "coarse-graining" procedure, the VBS 
states are transformed to a matrix-product form with 2- 
site block spins, characterized by the block-independent 
matrices Since the 2-site block spin states 
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|Gi2,M{|) form a complete orthogonal set, the null space 
in the reduced density matrix of a 2-site block is given 
by those IR channels with i?[<5i2'A^S] = Q. 

The null space for more than two adjacent sites is no 
longer easily visualized. However, the blocking process 
of g matrices can be proceeded without any fundamental 
difficulties (See Fig. El). In Sec. Ill B, we will study 
the spin-1 fermionic VBS state by using this powerful 
method. 

The uniqueness of the VBS ground states of the con- 
structed Hamiltonians has to be further clarified. Gen- 
erally, the ground state degeneracy will occur if there 
exists another state with a larger null space in the re- 
duced density matrix of the present interaction range. 
To lift the degeneracy, one can locate the null space in 
an extended range by blocking more spins and modify 
the Hamiltonian correspondingly. To justify the unique- 
ness, it is helpful to numerically diagonalize an open 
chain Hamiltonian with several lattice sites. If the nu- 
merically calculated ground-state degeneracy is D^, i.e., 
the ground states are all contained in the matrix-product 
form, one can prove the uniqueness of VBS ground states 
by a mathematical induction method. The basic idea is 
to assume that the VBS states I^'q./?) in Eq. (O are the 
only ground states of a projector Hamiltonian H{N) of 
an open chain with N sites. Then, the ground states of 
an open chain with TV -I- 1 lattice sites can be written as 
the following superposition of I^'q,/?) and \m]y^i): 

l*iV+l> = ^/3",m«+i |*a,/3) «) |tojv+1> . (11) 

Q/S.mjv + i 

The vectors Itoat+i) on the site A^-|- 1 decouple from the 
excited states of H{N) because such a coupling do not 
gain energy from the projector Hamiltonian H{N + 1). 

Now we can change the notation Wpa.mN+i = 1^/3™"^^' 
and then | ^E* n+i) can be immediately written in a matrix- 
product form 

X |mi • ■■ITlN+l) , (12) 

where D x D matrix can be fully determined 

because I^'at+i) are the zero-energy ground states of 
H{N + 1). After solving this eigenvalue problem, one 
can find that I^'at+i) can be written in the form of Eq. 
(ID). This final step completes the mathematical induc- 
tion proof. For periodic boundary conditions, the VBS 
ground state should be a linear combination of the 
states in Eq. ^ and be annihilated by the extra projec- 
tors acting on the two ends of the chain. This will lead 
to the VBS ground state in the form of Eq. 

Actually, those matrix product states with lin- 
early independent I^'q./s) in a finite open chain satisfy the 
so-called injective property i^^i"^^'^^ The injectivity of the 
matrix product states not only ensures the ground-state 



uniqueness of the parent Hamiltonian, but also guaran- 
tees the existence of an energy gap and the exponentially 
decaying correlation functions of local operators. 



III. GENERAL VBS STATES FOR QUANTUM 
INTEGER SPIN CHAINS 

In Sec. II, we set up a framework to study the VBS 
states with a built-in Lie group G. To test these abstract 
formalism, we begin with the quantum integer spin-S* 
chains and consider two classes of VBS states. In these 
two VBS classes, the virtual particles transform under 
spin-J representations and 2'^-dimensional S0{2S + 1) 
spinor representations, respectively. Toward the first 
class, Sanz et. al^ have explored 5'C/(2)-invariant two- 
body Hamiltonians which have such states as their eigen- 
states. In the present work, we further require that 
these VBS states are unique ground states of the par- 
ent Hamiltonians. However, the price we usually have 
to pay is to include multi-spin interactions in the par- 
ent Hamiltonians. This situation will be treated for the 
spin-1 fermionic VBS state in Sec. Ill B. For the second 
class, we will show that these states are equivalent to the 
S0{2S + 1) symmetric matrix product states introduced 
in our previous worki^ However, the present formalism 
explains the origin of the emergent SO{2S symmetry 
in these VBS states and shows that their edge states are 
50(25* -t- 1) spinors. Therefore, the two VBS classes are 
sharply distinct from each other for 5 > 3. 



A. Spin-5 VBS states with virtual spin-J particles 

As a warm up, let us apply the formalism in Sec. II to 
the spin-5 VBS states with two virtual spin-J particles 
in each site. It is well-known that the product of two 
spin-J representation 

J(g)J = 0©l®---®2J (13) 

always contains a singlet and the physical spin-5 repre- 
sentation if J > 5/2. After replacing the SU{2) Clebsch- 

Gordan coefhcient P^"^' = (5, to| J, a; J, /?) in Eq. ([T|) and 
the spin-J valence-bond singlet 

a=-J 

in Eq. the (2 J + 1) x (2 J + 1) matrix ^I'"' in Eq. 
(HI) can be written as 

A^"-\ ^Y.{-iy+0{S,m\J,a;J,-f3)\J,a) {J,(3\, (15) 

where — J < a,p < J. These A^"^^ matrices are rank 
5 irreducible spherical tensors and satisfy the following 
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commutation relations: 



B. Spin-1 fermionic VBS states 



[J±,A[™1] = v/(^Tm)(5±m + l)^['"±il, (16) 

where J± and Jz generate the spin- J representation of 
the SU{2) algebra, 



J± = V(J T ± a + 1) I J, a ± 1) (J, a\ , 

a 

J, = ^a|J,a)(J,a|. (17) 



For the celebrated VBS states of the AKLT models, 
i.e., the case of virtual spin J = S/2, we can also use 
the Schwinger boson representation to express the VBS 
states. In the Schwinger boson language, the spin oper- 
ators are expressed by 

S+ = alb,, S- - bla,, St = {a\a, - h%)/2, (18) 

with a local constraint a\ai+b\bi ~ 2S. Then, the integer 
spin-S' VBS states of the AKLT models in a periodic 
chain are expressed asii 



|AKLT)=n(«Ifcl+i-&IaI+i) 



(19) 



where |v) is the vacuum with no particle occupation. The 
matrix product form of these VBS states are obtained 
by Totsuka and SuzukiSi. Since the null space of two 
neighboring sites is the total bond spin S + 1,...,25' 
channels, the VBS states in Eq. are exact ground 

states of AKLT Hamiltoniansi^iii 

2S 

i?AKLT = 5I Yl JsTPsAhl + l). (20) 
i St=S+1 

where all > and Ps^ {h j) is the projection operator 
on total bond spin channel St- These SU{2) invariant 
projection operators can be written as polynomials of 
spin-exchange interactions • Sj up to 2S powers 



2S 



psAij)^ n 



2S,-Sj+2S{S+l)-S'{S' + 1) 



s'=o, 

S'^St 



SriSr + l)- S'iS' + 1) 



(21) 

For S/2 < J < S cases, at first glance, the null space 
of two neighboring sites is given by the total bond spin 
channels 2 J -I- 1, ... , 25", which is smaller than the VBS 
states of AKLT model. According to Sec. II B, one 
may conclude that next-nearest neighbor interactions are 
needed to be construct their parent Hamiltonians. How- 
ever, there is an exception: the VBS states with S — 2 
and J — 3/2. We will discuss this special case in Sec. Ill 
C. 



Now we consider S ~ 1 VBS state with virtual spin 
J — 1, which belongs to the class in Sec. Ill A. Since 
S" = 1 is the only antisymmetric product of two virtual 
J — 1 particles, one can use the fermionic statistics to 
implement the projection onto the physical S* = 1 sub- 
space. Thus, the physical S — 1 states are written as 



|l)^44|v), |0)=cIcLjv>, \-l) = cy_,\v). (22) 

The SUi2) spin operators are Sf = J2l,u=i ^l^S^a^ 
(a = X, y, z), where 5° are the usual 3x3 spin-1 matrices. 
The total spin = 2 on each site is imposed by a local 
constraint Yfn=i 4/^^^^ = 2. 

In terms of these fermionic variables, the S ^ 1 VBS 
state with virtual spin J = 1 can be exactly written as 

TV 

l*i> = n(4i4+i,-i-4o4+i,o+4,-i4+i,i) 1^) ' (23) 

which has a matrix product form in Eq. Q with 



t J 



^[1] 




MM 




(24) 



Following the method in Sec. II B, we can construct 
the parent Hamiltonian for this fermionic VBS state. To 
locate the null space, it is sufficient to block three suc- 
cessive spins. The tensor decomposition of three spin-1 
representation yields 

l(8)l(g)l = (0©1©2)«)1 

= 1©0®1'©2©1"®2'®3, (25) 

which provides a natural basis for block spins. In this 
basis, the block states can be denoted by \Si2;S,M), 
where S and M are total spin and magnetic quantum 
number of the three sites, S12 is the total spin of the 
first two sites. For the representations 1' and 2, we have 
512 = 1- For the representations 1" and 2', 5*12 = 2. 
For the representations and 3, the index S'12 can be 
suppressed and does not lead to misunderstanding. 
After blocking the three spins, wc obtain 



515253 = 



E 

mim2m3 



S,M S12 

[S.M] 



(26) 



where the 3x3 matrices Cg^^ are given by 
(S'i2,mi -f m2|l,mi; 1,TO2) 



C 



[S,M] 
S12 



X (S',M|S'i2,mi +m2;l,m3)^['"iU["^U["^l. 
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The matrices C^Si^^ '-^'^ calculated by using Eq. p4)) . 
Firstly, we find that C[°'°l ^ and Cl^-*^! = 0. This 
means that the spin-0 singlet state is contained but the 
spin-3 states are absent in every three-site block. The 



other nonvanishing matrices Cg_^ 
equations: 



[SM] 



satisfy the following 



C. 



C. 



[2,M] 



[2,M] 



:i,M] 



5 ^[13/] 



V5 [i,Af] 

—^0 ■ 



(27) 



According to Eq. , the unnormalized states contained 
in the 3-site block 5i(?233 are one spin-0 state, three spin- 
1 states 

4 |0; 1, M) - V3 |1; l,M) + V5 |2; 1, M) 

with — 1 < M < 1, and five spin-2 states 

|1;2,M) + V3|2;2,M) 

with — 2 < M < 2. By using the Gram-Schmidt orthogo- 
nalization method, we find that seven spin-3 states |3, M) 
with —3 < M < 3, five spin-2 states 

|02.Af) = V3|1;2,M)-|2;2,M) 

with —2<M < 2, six spin-1 states 

\/3|0;l,M)+4|l;l,M), 



91, m; 
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|^i,m) = 4|0;l,Af)-V3|l;l,M)- — |2;1,M) 

with — 1 < Af < 1 span the null space in the reduced den- 
sity matrix of the 3-site block. Therefore, the 3-site pro- 
jector Hamiltonian for which the spin-1 fermionic VBS 
state is the zero energy ground state is thus given by 

/i = A3 5] |3,A/)(3,M| + A2 l^2,Af) (</'2,m| 

|A/|<3 |M|<2 

+ ^ (Ai |(/.i,Af) (01,Af| + A'l |(pi,A/) ((pi,A/|)(28) 
|Af|<l 

where all A3, A2, Ai, A'^ > 0. 

It is interesting to compare the spin-1 fermionic VBS 
state with the spin-1 bosonic VBS state of AKLT model. 
In the fermionic VBS state, the two-point spin corre- 
lation function decays exponentially with a correlation 
length ^ = l/ln2, longer than that for the AKLT model 
= l/ln3). Besides the obvious difference of the edge 
spin representation in an open chain, we can also see a 
sharp difference by computing the non-local string order 
parameter- 
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0{e)= hm (5,fnexp(«^^.^)5| 



(29) 



By using the transfer matrix method^^, we arrive at the 
result 0{9) = ^sivP'O for the fermionic VBS state. For 



CD 
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FIG. 3: The non-local string order parameter 0{6) as a func- 
tion of the spin- twist angle 6 for the fermionic VBS state (solid 
line) and the bosonic VBS state (dashed line) of the AKLT 
model. 



comparison, the values of the non-local order parameter 
0{9) for both the spin-1 fermionic VBS state and bosonic 
VBS state of the AKLT model are plotted in Fig. © . For 
the spin-1 bosonic VBS state, O{0) reaches its maximum 
at = TT, which is reduced to the den Nijs-Rommelse 
string order parameter characterizing the hidden antifer- 
romagnctic order in the AKLT VBS state. However, in 
the fermionic VBS state, 0{9) has a minimum ©(tt) = 
and the maximum value 1/9 at both 9 = tt/2 and 37r/2. 
This signifies that the hidden antiferromagnetic picture 
totally breaks down in the S = 1 fermionic VBS state. 
How to describe such a state has not been clear so far. 



C. VBS states with an emergent 50(25' + 1) 
symmetry 

In Sec. Ill A, we have mentioned an exceptional exam- 
ple: S = 2 VBS states with virtual spin J = 3/2. Besides 
the absent bond total spin S't = 4 channel of neighbor- 
ing sites, there is a new forbidden channel St = 2 in this 
VBS statcj^ Therefore, its parent two-body Hamiltonian 
can be written as 



H 



[JlP2(i,«+l) + J2P4(i,i+l)] 



(30) 



with Ji, J2 > 0. According to Eq. (PTjl . the projector 
Hamiltonian (|5D|) can be rewritten as 



H 



E 



3 J2 - 8OJ1 



+ - 



84 

10 Ji + J2 
60 



S,- • S-, 



9 J2 - 40 Ji 
360 



{Si ■ Sj 



(31) 



In fact, the 5 = 2 VBS state with virtual spin J = 3/2 
has a hidden 5*0(5) symmetry and its matrix-product 
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form was studied by Scalapino et al^ in an 5*0(5) sym- 
metric ladder system of interacting electrons. 

The quantum spin-2 Hamiltonian in Eq. (|30|) be- 
longs to a new class of exactly solvable quantum integer- 
spin chains introduced by the first three of us very 
recently)^ The ground states of these Hamiltonians are 
SO{2S 1) symmetric matrix-product states and exhibit 
hidden topological order. For 5 = 1, the 50(3) symmet- 
ric matrix product state becomes the VBS state of spin-1 
AKLT model. For 5 = 2, the 50(5) symmetric matrix 
product state is the 5 = 2 VBS state with virtual spin 
J = 3/2. However, it was not clear whether this fam- 
ily of matrix product states has a valence-bond picture 
for 5 > 3. By using the framework in Sec. II, we will 
show that there is indeed a VBS picture for these matrix 
product states. Actually, the virtual particles in these 
VBS states transform under the 2'^-dimensional spinor 
representation of SO{2S + 1). 

It is convenient to promote the symmetry of the sys- 
tem and demand the spin-5 states on each site transform 
under the (25 -f l)-dimensional vector representation of 
50(25 + 1). The tensor product of two 50(25 + 1) vec- 
tors on the adjacent sites can be decomposed as 

25 + 1 ^ 25 + 1 = 1 @ 5(25 -I- 1) ® 5(25 + 3) , (32) 

where the number above each underline is the dimen- 
sion of the corresponding IR. These 50(25 -I- 1) IRs can 
be directly related to SU{2) integer-spin IRs. Here 1 
is the symmetric spin singlet, while the antisymmetric 
channel 5(25 + 1) and the symmetric channel 5(25 + 3) 
correspond to the total bond spin 5t = 1, 3, . . . , 25 — 1 
and 5t = 2, 4, . . . , 25 states, respectively. Therefore, the 
SO{2S + 1) bond projection operators can be expressed 
using the spin projection operators as 

s 

Ps{2S+i) {i,j) = ^PsT=2i-i{i,j), (33) 
1=1 
s 

Ps(2S+3) {hj) = Y.PSr=2l{hj)- (34) 
1=1 

On each lattice site, the SO{2S + 1) vectors can 
be formed by tensor decomposition of two virtual 2"^- 
dimensional spinors 

2^«2£^e(,/,,), (35) 

q—O — — 

where (25^1) = q\(2S-q+iy. • Note that q — O and g = 1 in 
Eq. ((35|l correspond to singlet representation and (25 + 
l)-dimensional vector representation, respectively. For 
5 = 1, Eq. ([55)1 recovers the well-known decomposition 
2(g)2 = le3of two spin-1/2 spinors. For 5 = 2, Eq. 
(|35|) can be interpreted as the decomposition 4 ® 4 = 
1 ® 5 ©10, where the 50(5) spinors can be viewed as 
spin-3/2 variables because 50(5) ~ 5p(4). However, the 



50(25 + 1) spinors in Eq. JSSl) for 5 > 3 do not have 
SU{2) spin counterparts. 

Following the discussions in Sec. II, the 50(25 + 1) 
symmetric VBS states can be constructed by combining 
the virtual spinors on the neighboring sites into valence- 
bond singlets. By comparing Eq. ([311) and Eq. (|35p . 
one finds that the IR channel 5(25 -I- 3) for any two 
neighboring sites is absent in these VBS states. Here 
an interesting observation is that those IR channels with 
q > 2 in Eq. ([55]) are actually absent for two adjacent 
sites due to the projection of two virtual spinors onto 
the physical vector representation in each site. There- 
fore, the SO{2S + l)-invariant parent Hamiltonian for 
the SO{2S + 1) symmetric VBS states is given by 

H = Y, Psi2S+3) {'.i + l)- (36) 

i 

Since the null space of these VBS states is the non-zero 
even total spin channels, we can extend the SO{2S + 
l)-invariant parent Hamiltonian to the following SU{2)- 
invariant quantum integer-spin Hamiltonian 

s 

H = Y,Y.Ji^ST=2i{^,i + l), (37) 

i 1=1 

with aU J; > 0. 

Actually, the SO{2S + 1) symmetric VBS states are 
equivalent to the matrix product states studied in Refi^. 
In the present VBS form, the origin of emergent 50(25-|- 
1) symmetry and the 2^^ edge states on each boundary of 
an open chain are quite clear. Although the edge degrees 
of freedom in 5 = 1 and 5 = 2 cases can be viewed as 
SU{2) spin variables, they transform under SO{2S + 1) 
spinor representation for 5 > 3 cases. It is interesting to 
compare these SO{2S + 1) symmetric VBS states to the 
spin-5 VBS states formed by virtual spin J = {2^ — l)/2 
in Sec. HI A. Although they are both unique in a pe- 
riodic chain and 4'^ fold degenerate in an open chain, 
their distinct edge states show that they belong to two 
different topological classes. These explicit examples im- 
ply that the ground state degeneracy is not sufficient to 
characterize the topological ordered states. 

IV. 50(5) SYMMETRIC VBS STATES 

So far, we are restricted to the case of SU{2) integer- 
spin in each site. Actually, the method discussed in Sec. 
II can be applied for a general Lie group G, we thus 
move on to 50(5) Lie group, where the physical states 
transform under 50(5) IRs. 

The 50(5) Lie algebra has 10 generators i"'' (1 < a < 
6 < 5), satisfying the commutation relations 

L'^] = liSadL"'' + SbcL'"' - SacL""" - SmL^- (38) 

Mathematically, the IRs of 50(5) are labeled by two inte- 
gers (p, q), with p > q > 0. For the (p, q) representation 
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TABLE I: Several irreducible representations of the SO{5) Lie 
group. 



Representation 



Dimension 



Casimir charge 



(0,0) 
(1-0) 

(1,1) 

(2,0) 
(2,2) 
(3,1) 
(4,0) 



1 
4 
5 
10 

14 
35 
35 





5/2 
4 
6 

10 
12 

16 



of 50(5) Lie group, the dimensionality d{p, q) and the 
Casimir charge C {p, q) are given by^ 

d{p,q) = (l + 5)(l+p-g)(l + |)(l + ^),(39) 



C{p,q) = Y.iL^'f = 



9 ? 

p^ q'^ 



2p + g, 



(40) 



a<b 



respectively. The dimensionality and Casimir charge for 
the simplest 50(5) irreducible representations are listed 
in Tab. m 



A. Bosonic 5*0(5) VBS states 

We begin with the 10-dimensional (2, 0) adjoint rep- 
resentation of 50(5). The bosonic 50(5)/5p(4) VBS 
state of this system was first considered by Schuricht and 
Rachel)^ Their strategy is to construct the (2, 0) adjoint 
representation by two virtual particles transforming un- 
der the (1,0) spinor representation, 



(1,0)® (1,0) = (0,0) ©(1,1)® (2,0). 



where (0, 0) and (1, 1) are antisymmetric and (2, 0) is the 
only symmetric product representation. Therefore, one 
can obtain the physical (2,0) adjoint representation by 
endowing bosonic statistics to the virtual (1, 0) spinor 
particles. This is analogous to the SU{2) Schwinger 
boson representation which symmetrizes two spin-1/2 
spinors to construct a spin-1 representation. Using the 
4-component 50(5) Schwinger bosons, the 50(5) gener- 
ators in Eq. can be defined by 



1 4 

L"-'' = - b'^T°-''b 



(42) 



where F"^ = [T'',T^]/2i and 



t1,2,3 





—id 



/ 
/ 



/ 

-/ 



(43) 

For the (2, 0) adjoint representation with ^Jt^M = 2, 

the 10 states in a bosonic language are shown in the (2, 0) 
weight diagram in Fig. H) After a rotation by 45°, this 



h\h\ |v) 



-1 



r\2 
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L34 

hXhX |v) 



,t^2 



b\h\ |v) 



h\bl |v) 



blh\\ 



b\h\\ 



FIG. 4: Weight diagram and the bosonic realization of the 
(2,0) adjoint representation of SO{5). There is a two-fold 
degeneracy with L^^ = L'^'^ = 0. 



weight diagram is identical to that given by Schuricht and 
Rachel. Here we choose the Clifford algebra generated 
by the F matrices to define the 50(5) generators. The 
advantage of our convention is to find an interesting non- 
local hidden string order in the (2,0) bosonic VBS state 
below. 

The (2, 0) bosonic 50(5) VBS state is formed by con- 
tracting two (1,0) spinors on neighboring sites into a 
valence-bond 50(5) singlet. Its wave function can be 
written compactly as 



i*2)-n(E^L^A..^i+i.jiv), 



(42) where the antisymmetric matrix R is given by 



R 







-iay 
with the following properties: 



i?2 = -1, i?t = = i?^ = -R, 



RT'^R-' = (F°)' , Rr''''R-'^ = -(F"") 



2b\T 



(44) 



(45) 



(46) 



The tensor product decomposition of two neighboring 
(2, 0) adjoint representations is written as 

(2, 0)® (2, 0) = (0, 0)® (1, 1)® (2, 0)® (2, 2)® (3, 1)® (4, 0). 

(47) 

In the (2, 0) bosonic VBS state, a valence-bond singlet 
of two virtual (1,0) spinors are created and therefore 
the two adjacent sites can only transform (0,0), (1,1), 
and (2,0) representations according to Eq. ([1T|) . Con- 
sequently, \'^2) is an exact ground state of the projector 
Hamiltonian 

^ = II ['^1^(2,2) («,« + !) + J2P(3,1) (*,«+!) 

i 

+^3^(4,0) + 1)] , (48) 
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where Ji, J2, J3 > and P(2,2),^(3,i),^(4,o) are projec- 
tors onto the (2, 2), (3, 1), (4, 0) representations, respec- 
tively. 

Furthermore, the (2, 0) bosonic VBS state contain a 
well-defined hidden string order. This can be observed 
in its matrix-product wave function with the local matrix 



/ |0,0) V2|l,l) 

-x/2|-l,-l) -|0,0) 
|0,-1) |1,0) 

V -1-1,0) -|0,1) 



|0,1) |1,0) \ 

-1-1,0) -|0,-1) 
|0,0)' V2|l,-1 
-V2|-l,l) -|0,0)'/ 



where we take |0,0) = b\bl |v) and |0,0)' = blbl |v). In 
both of the toi and m2 channels, it can be shown that 
|m,,) (?7 — 1,2) has a hidden antiferromagnetic order. In 
other word, the states of m,, — 1 and —1 will alternate in 
space if all the — states between them are ignored. 
A typical configuration of this state is given by 

mi: •••OTOOiTiOOOtOiOt--- 
m2: •••OTOOOiTOiOOTiOO--- 

where (T,0,i) represent \m) = (|1), |0), |-1)). This 
hidden antiferromagnetic order is in analogy with the 
spin-1 VBS state of AKLT modcl^^ and its 50(25 + 
1) generalization'^^. To characterize this hidden anti- 
ferromagnetic order, one can generalize the den Nijs- 
Rommelse string order parameters as 



lim 



(if W exp{iTTLf)L 



(49) 



In fact, the non-local string order parameters for the Car- 
tan generators introduced by Schuricht and Rachel'^^ are 
combinations of our O^^ and O^^. The advantage of 
our convention is that the string order parameters in Eq. 
(I49p clearly reflect a hidden antiferromagnetic order in 
the (2,0) bosonic SO{5) VBS state. The value of these 
string order parameter can be obtained by a probabil- 
ity argument. These non-local string order parameters 
should all be equal to each other because the VBS state 
preserves 50(5) symmetry. Thus, we only need to eval- 
uate the value of O^^ by considering mi channel. The 
role of the non-local string phase factor in Eq. (I49p is to 
correlate the finite spin polarized states in the mi chan- 
nel at the two ends of the string. If nonzero mi takes 
the same value at the two ends, then the phase factor is 
equal to 1. On the other hand, if a nonzero mi takes 
two different values at the two ends, then the phase fac- 
tor is equal to —1. Thus, the value of O^^ = 9/25 is a 
square of the probability of the nonzero mi = ±1 ap- 
pearing at the ends of the string. Correspondingly, a 
generalized Kennedy- Tasaki unitary transformation can 
be designed according to Ref.^'' and we expect that the 
50(5) symmetry of the original Hamiltonian is reduced 
to {Z2 X ^2)^ under such a non-local transformation. The 
non-local string order parameters in Eq. (|49p for the Car- 
tan generators will be transformed to two-point correla- 
tion functions, which properly characterize the hidden 



{Z2 X ^2)^ symmetry breaking. Thus, in this state, the 
non-local string order and the 16-fold degeneracy in an 
open chain can be viewed as natural consequences of a 
hidden (Z2 x ^2)^ symmetry breaking. 

In fact, the bosonic 50(5) VBS state of (2,0) adjoint 
representation can be generalized to the totally symmet- 
ric (p, 0) representation with even p. Namely, {p, 0) rep- 
resentation for even p can be constructed by two (p/2, 0) 
representations. However, there is an alternative way to 
take the advantage of a generalized Schwinger boson rep- 
resentation. Using the generalized Schwinger boson rep- 
resentation, the {p, 0) representation in each site can be 
constructed by symmetrization of p spinors and the local 
constraint is now replaced with ^Ii^ij. ~ P- Thus, 

the bosonic 50(5) VBS states for {p, 0) representation 
can be written as 



(50) 



In an open chain, there are fractionalized edge states 
transforming under (p/2, 0) representation. Once a CP^ 
coherent state representation^^ is used, the (p, 0) VBS 
states have the Jastrow form. They are analogous to the 
fractional quantum Hall states in CP^ spaced at filling 
fraction v = 2/p, in the same sense as the resemblanceii 
between VBS states of AKLT model and the fractional 
quantum Hall states in spherical geometry. 

Since the tensor product decomposition of two (p, 0) 
representations is given by 



(p,0)®(p,0) = j]^(fc + /,fc-/). 



(51) 



fc=0 (=0 



and p/2 valence-bond singlets are created between ad- 
jacent sites in l^s), the only finite projections on two 
adjacent sites are given by 



p/2 k 



{p/2, 0) ® {p/2, 0) = Y.Y.^^ + ^' ^ " 0- 



(52) 



fc=o 1=0 



Thus, the null space of the 2-site reduced density ma- 
trix are given by a sum of the representations written as 

X]fe=p/2+i X]f=o(^"'"'' fc — and the corresponding parent 
Hamiltonian of j^Ps) is given by 

p k 

,k-l)P{k+l,k-l){^,^ + '^)^ (53) 

i +1 /=o 

where all J{k+i,k^i) > and P(k+Lk-i) is the projection 
operator onto the {k + l,k — I) representation states. 

With the help of the Casimir charge in Eq. the 
50(5) projectors can be written as polynomial functions 
of 50(5) generators. According to Eq. (|5ip . these pro- 
jectors satisfy a completeness relation 



^^P{k+i.k-i){h3) = 1- 



(54) 



k=0 1=0 
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Considering the two-site Casimir charge X]a<b(^i'^ + 
L'f'Y ^ we can write the S'0(5) Heisenberg interaction as 



Y^LfLf = i^^[C(fe + ;,fc-0-(p'+4p)] 

a<b k=0 1=0 

>^P(k+i,k-i){i,j)- (55) 
Using the properties of the projectors, we have 

p k 

2 



a<b 



k=0 1=0 



{k+l,k-l) 



Together with the completeness relation (|54|) . this for- 
mula can be inverted, so that each projector can be rep- 
resented by a polynomial function of 5*0(5) Heisenberg 
interaction Ea<h^f • 



L34 



i(ci + i4)(4-ici)|v> 



c\ + ^4)4 |v) 



(56) i(4 + *4)(4 + *4)|v) 



i(4_,4)(4_,4)iv) 



44 |v) ^(4-*4)4lv) 



44 |v) 



i(4-i4)(4 + i4)iv) 
i=(4-i-z4)4iv) 
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FIG. 5: Weight diagram and the fermionic realization of the 
(2, 0) adjoint representation of SO{^) Lie algebra. 



B. Fermionic SO{^) VBS state 

In this subsection, we present another way to construct 
the (2,0) adjoint representation, i.e., by using two (1, 1) 
vector representations, 

(1,1)® (1,1) -(0,0) ©(2, 0)0(2,2), (57) 

where the (2, 0) adjoint representation is antisymmetric 
and (0,0), (2, 2) are symmetric. This is because the or- 
thogonal groups have a general property that the adjoint 
representation is the only resulting antisymmetric chan- 
nel of two vector representations.^"* The simplest realiza- 
tion of this property is the 50(3) spin-1 case discussed 
in Sec. Ill B, where the antisymmetrization of two vec- 
tor spin-1 representations only yields the spin-1 adjoint 
representation. 

If we use the fermionic statistics to implement the an- 
tisymmetrization, the 10 states in the adjoint represen- 
tation can be written as 44 1^)' where 1 < a,6 < 5. 
Moreover, the 50(5) generators are defined by 

L'^' ^i{cicb-clca), (58) 

and a double occupancy constraint X]o=i — 2 can 
guarantee the adjoint representation in each lattice site. 
Using these fermionic variables, the (2, 0) weight diagram 
is shown in Fig. O 

Using the fermion variables, the (2,0) fermionic VBS 
state with two virtual (1, 1) vector 50(5) representations 
can be written as 

l*4)-n(E4Al,a)|v). (59) 

i a 

In an open chain, the edge spins transform under (1, 1) 
vector 50(5) representation, different from the (1, 0) 
spinor 50(5) representation in the (2, 0) bosonic 50(5) 



VBS state. Another interesting observation is that 
the perfect non-local string order presence in the (2, 0) 
bosonic 50(5) VBS state vanishes in the fermionic VBS 
state, because the string order parameter (|49)) for this 
state is found to be zero. In this sense, the bosonic and 
fermionic (2,0) VBS states can be viewed as 5*0(5) gen- 
eralizations of spin-1 VBS states of AKLT model and 
fermionic VBS state in Sec. Ill B. 

Finally, using two-body interactions, one can construct 
the parent Hamiltonian for this fermionic 5*0(5) VBS 
state. Since any two adjacent sites can only transform 
under (0,0), (2,0), and (2,2) representations, |5'4) is an 
exact zero-energy ground state of the projector Hamilto- 
nian 

H = ^ [ifiP(ia)(i,z + l) + if2P(3,i) («,* + !) 

i 

+K3P(i^o){hi + l)], (60) 

for Ki, K2, K3 > 0. The possible hidden order is still 
under investigation. 

V. CONCLUSION 

In conclusion, we have presented a general method to 
construct one-dimensional VBS states embedded with 
Lie group G and their parent Hamiltonians. This pro- 
vides examples that the topologically ordered states can 
be systematically generated in one dimension and are 
characterized by their edge states representations as well 
as their ground state degeneracy. 

For quantum integer spin-5' chains, there exists two 
topologically distinct families: (i) the virtual particles 
transform under SU(2) spin-J representations and (ii) 
the virtual particles are 5*0(25 + 1) spinors. In the first 
class, a new spin-1 fermionic VBS state is constructed as 
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an explicit example. Compared to the celebrated S = 1 
valence bond solid state of AKLT model, the fermionic 
valence bond solid state shows drastic differences on the 
edge states and hidden string order. For the second class, 
it has been shown that these valence bond solid states 
with an emergent SO{2S+l) symmetry are equivalent to 
the previously proposed SO{2S +1) symmetric matrix- 
product statesi^ The present formalism explicitly dis- 
plays that the edge states of an open chain transform 
under the SO{2S + 1) 2"^-dimcnsional spinor representa- 
tion. 

To generalize the VBS states in SU{2) symmetric 
quantum integer-spin chains, two types of VBS states 
with the SO{5) symmetry are considered, including (i) 
bosonic SO{5) VBS states formed by a symmetrization 
of two spinor representations in each site and (ii) a 
fermionic 5*0(5) VBS state with (2, 0) adjoint represen- 
tation formed by antisymmctrization of two vector rep- 
resentations. 



It can be expected that the ideas and formalism devel- 
oped in this work are very useful and can be generalized 
to the tensor product states (projected entangled pair 
states) for higher dimensional correlated systems;^ The 
understanding of the physical properties of these states is 
the first step to characterize higher dimensional topologi- 
cal states, which certainly deserves further investigations. 
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